Putting a Spin on the Gaffnian 



Simon C. Davenport, 1 Eddy Ardonne, 2,3 Nicolas Rcgnault, 4,5 and Steven H. Simon 6 

1 Rudolf Peierls Centre for Theoretical Physics, Oxford University, 

1 Keble Road, Oxford, 0X1 SNP, United Kingdom 
2 Nordita, Royal Institute of Technology and Stockholm University, 
Roslagstullsbacken 23, SE-106 91 Stockholm, Sweden 
3 Department of Physics, Stockholm University, AlbaNova University Center, SE-106 91 Stockholm, Sweden 
4 Department of Physics, Princeton University, Princeton, NJ 08544, USA 
5 Laboratoire Pierre Aigrain, ENS and CNRS, 24 rue Lhomond, F-75005 Paris, France 
6 Theoretical Physics, Oxford University, 1 Keble Road, Oxford, 0X1 3NP, United Kingdom 

(Dated: October 31, 2012) 

We characterize in detail a wave function conceivable in fractional quantum Hall systems where a 
spin or equivalent degree of freedom is present. This wave function combines the properties of two 
previously proposed quantum Hall wave functions, namely the non-abelian spin-singlet state and 
the non-unitary gaffnian wave function. This is a spin-singlet generalization of the spin-polarized 
gaffnian, which we call the 'spin-singlet gaffnian' (SSG). In this paper we present evidence demon- 
strating that the SSG corresponds to the ground state of a certain local Hamiltonian, which we 
explicitly construct, and, further, we provide a relatively simple analytic expression for the unique 
ground state wave functions, which we define as the zero energy eigenstates of that local Hamilto- 
nian. In addition, we have determined a certain non-unitary, rational conformal field theory which 
provides an underlying description of the SSG and we thus conclude that the SSG is ungapped in 
the thermodynamic limit. In order to verify our construction, we implement two recently proposed 
techniques for the analysis of fractional quantum Hall trial states: The 'spin dressed squeezing 
algorithm'; and the 'generalized Pauli principle'. 



I. INTRODUCTION 

The study of the fractional quantum Hall effect 
(FQHE) of electrons in semi-conductor materials^ and, 
potentially, of bosons in rotating Bose gasses^ or in arti- 
ficial gauge fields 3 -, has been a fertile ground for the con- 
ception of topological phases of matter, which exhibit ex- 
otic and exciting theoretical properties^ In this paper we 
shall discuss the construction of a FQHE wave function 
exhibiting a 'multi-component' degree of freedom such 
as a spin 5 -, valley^, iso-spin£, layer— or sub-ban d 10 ' 11 . 
The wave function describes the low-energy behaviour of 
a system possessing a vanishing energy gap to quasihole 
excitations in the thermodynamic limit, and is thus not 
itself a candidate to describe a topological phase. Never- 
theless the study of this wave function is an interesting 
exploration that enhances the understanding of the con- 
struction of 'multi-component' FQHE wave functions in 
general. 

An ideal theoretical description of a FQHE state com- 
prises at least the following three ingredients: A local 
Hamiltonian which describes the ground state and exci- 
tation spectrum of the Hilbert space; relatively simple an- 
alytic wave functions which are the highest density zero 
energy states corresponding to the Hilbert space of that 
Hamiltonian; and a (rational) 2D conformal field theory 
(CFT) which generates these wave functions. If a plasma 
analogy is available, then a gapped state is associated 
with the analogous plasma being in a screening phased 
In general, the quantum Hall Hilbert space is built from 
a basis of monomials in the complex particle coordinates 
z%. In certain instances — the Laughlin series, the Rcad- 



Rezayi series, the gaffnian and Haffnian wave functions, 
the Halperin wave functions and the non-abelian spin 
singlet (NASS) states^— — it has been possible to i) 
determine simple analytical expressions for the polyno- 
mial wave functions, and ii) construct a local Hamilto- 
nian whose zero-energy eigenstates are in one-to-one cor- 
respondence with those polynomial wave functions. A 
key feature in each of these special cases is that the wave 
functions are uniquely defined by their vanishing proper- 
ties. 

Changing wave functions describing a FQHE state in a 
seemingly innocuous way, may in fact have severe conse- 
quences. The wave function of the Moore-Read state 
for spin-less bosons vanishes quadratically when three 
of the constituent bosons are coincident. The gaffnian 
wave function is obtained by changing this behaviour, 
such that the wave function vanishes as a third power, 
when three constituent particles are coincident^ This 
simple change, however, results in a non-unitary, com- 
pressible wave function, which does not describe a topo- 
logical phase. Nevertheless non-unitary wave functions 
such as the gaffnian are still of interest as they are 
thought to correspond to critical points between other 
unitary, incompressible topological phases. This scenario 
is well understood in the case of the Haldane-Rezayi^ . 
wave function, which describes the phase transition be- 
tween a d-wave spin-singlet phase and a strongly paired 
stat e 15 ' 21 . Similar scenarios have been suggested for the 
gaffnianJ^ It is noteworthy that the gaffnian wave func- 
tion has large overlap with an incompressible composite 
fermion state thus suggesting that the gaffnian is a criti- 
cal point between the composite fermion phase and some 
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other phasei ' 22 — 

In this paper, we present evidence for a quantum 
Hall wave function, which inherits properties from both 
the NASS state — a unitary, spin-singlet quantum Hall 
state — and the non-unitary gaffnian. Such a wave func- 
tion has also been considered in Ref. [2ft Most notably, 
we argue that the wave function can be written as a spe- 
cial polynomial and that it also corresponds to the high- 
est density zero energy state of a certain local Hamil- 
tonian. We christen it the 'spin-singlet gaffnian' (SSG). 
(Throughout this work, we shall describe how the wave 
function corresponds to a 'spin' degree of freedom, how- 
ever note that all of what follows applies equally well to 
any other of the other type multi-component degree of 
freedom such as valley or layer index.) 



the proposal of a model local Hamiltonian. First, we 
shall explain how this local Hamiltonian is constructed 
and then we shall present two alternative methods to 
count the number of zero-energy modes of the proposed 
Hamiltonian. The first method employs the spin-dressed 
squeezing algorithm (which can also be used to obtain the 
wave functions) and second method employs the general- 
ized Pauli principle. In Sec. lIIII we shall present numerical 
evidence demonstrating that the zero-modes generated 
by diagonalizing the local Hamiltonian presented here are 
in precise agreement with the zero-modes generated by 
the squeezing algorithm and the counting obtained from 
the generalized Pauli principle. 

A. Pseudopotential Construction of Hamiltonian 



Statement of Results 

Our first key result is a model local Hamiltonian de- 
scribing the SSG wave function. In Sec III Al we shall 
explain how it can be written in terms of generalized Hal- 
dane pseudopotentialsj^ir— In Sees. Ill Bl and lll Cl we give 
an overview of how both the squeezing algorith m 29 ' 30 and 
the generalized Pauli principle for spinful states^ apply 
to the SSG wave function. In Sec. IIIII we shall present 
the results of numerical exact diagonalization of the SSG 
Hamiltonian. We have determined that the counting of 
zero energy eigenstates in the quasihole spectrum and 
also the entanglement spectru m 31 ! 32 matches the count- 
ing predicted by the squeezing algorithm. The counting 
of zero-energy eigenstates is also found to be identical to 
the result generated by the spinful version of the gener- 
alized Pauli principle. 

Our second key result is an analytic form for the 
ground state wave function. The SSG wave functions 
are constructed from conformal blocks in a CFT asso- 
ciated with the semi-direct product of non-unitary min- 
imal models expressible as M(3,5) x M(5,7) (see also 
Ref. HI). We shall discuss the CFT in section El In 
a bosonic incarnation, the SSG ground state occurs at 
filling factor v = 4/5 (and on the sphere, it has a shift 
6 = 3). The proposed analytic construction for the SSG 
is consistent with the vanishing properties required by 
the CFT and the ground state of the Hilbcrt space gen- 
erated by squeezing. We have determined that our pro- 
posed ground state wave function is the unique, highest 
density zero energy ground state of our model Hamilto- 
nian at v = 4/5 and <5 = 3. We shall discuss an explicit 
form of the ground state wave function in Sec. [Vj 



II. MODEL HAMILTONIAN AND THE 
CHARACTERIZATION OF ITS ZERO MODES 

In this section we shall describe in detail, three inde- 
pendent methods by which we are able to study the SSG 
wave function. The cornerstone of our argument lies with 



In prior investigations, most notably for the Laugh- 
lin and Moore-Read wave functions, it has been de- 
termined that trial quantum Hall wave functions cor- 
respond to unique, zero-energy ground states of certain 
model Hamiltoniansi 26 ' 33 Often, these model Hamilto- 
nians are most simply expressible in terms of Haldane 
pseudopotcntials^, and their generalization s 27 ! 28 . Given 
a system with a certain M-body interaction potential 
V(z\, ...Zm)i the action of a pseudopotential is to project 
out a particular component of that interaction. Those 
components are labelled by a convenient set of quantum 
numbers, which describe all possible few-particle inter- 
actions. The vector space of few-particle interactions 
is spanned by specifying both the relative angular mo- 
mentum L and, if we consider a spin degree of freedom 
as well, the spin quantum number S also. (Note that 
in general it is necessary to further distinguish between 
distinct interaction components with the same L, S. In 
other words, there exists in general a sub-vector space 
of dimension <Il,s for each L and S sectori 28 ' 34 In this 
paper, however, we shall only be concerned with sub- 
spaces of dimension d,L,s=l or 0, and so, for clarity, we 
omit any additional notation.) Pseudopotcntials spec- 
ify projection operators in the Hamiltonian and are thus 
labelled by two distinct sets of such quantum numbers. 
The pseudopotential V^ S . L , S , is expressed in terms of the 
vector space \L,S) and the M-body interaction potential 
V(zi, ...z M ) as 

V™ S . L , S , = {L,S\V(z u ...z M )\L',S') . (1) 

The general form of the Hamiltonian is then given by 

H= \L,S)V L % L , S ,(L\S>\. (2) 

L,S.L'S' 

Such a basis of pseudopotcntials is particularly con- 
venient when the interaction potential is rotationally 
and/or spin rotationally invariant, since in that case the 
pseudopotentials are diagonal in L, S, or both. 

At our convenience, we can pick certain special many- 
body interactions (such as delta function type interac- 
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tions) for which only a small set of pscudopotcntials re- 
main non-zero. The impact of specifying a positive value 
of a given pscudopotcntial in a model Hamiltonian is to 
assign energy to the corresponding component of the in- 
teraction, therefore, if such a component is present in a 
given trial wave function, then that wave function will 
not be a zero energy state of our Hamiltonian. Con- 
versely, if a component is not present in a given trial 
wave function, then we can include the corresponding 
pscudopotcntial in the Hamiltonian without introducing 
an extra energy cost. In this way we can tailor the Hamil- 
tonian to correspond to the desired properties of a given 
trial ground state wave function (these properties might 
come, for instance, from a CFT description of the state, 
see SecHVl). 

To give an example, briefly, it is known that the 
Moore-Read wave function for bosons at v = 1 is 
the unique, highest density zero energy ground state of 
a spin-polarized three-body contact interaction^ The 
space of pscudopotcntials for spin-polarized three-body 
interactions is spanned by a relative angular momentum 
L. S takes only its maximal value for a three-body in- 
teraction, S = 3/2i21 (Note a slight abuse of notation 
here: We are denoting bosons as if they are spin- 1/2 ob- 
jects, e.g. three bosons can have S = 3/2 and 5=1/2 
interaction channels. In actual fact we are describing sys- 
tems such as cold atomic gasses for which bosons can be 
engineered to have internal two-state degrees of freedom, 
and we have simply mapped that two-state degree of free- 
dom onto a spin degree of freedom.) In the language of 
pseudopotentials, the Moore- Read state corresponds to 
^o 3 3/2-o 3/2 being positive and all other pseudopotentials 
being zero and, in addition, the corresponding Hilbcrt 
space is restricted to only spin-polarized sectors. Along 
similar lines, the spin-polarised gaffnian wave function 
for bosons corresponds to V^ 3 3 y 2 . 3 ^ 2 and 3 / 2 . 2 3 / 2 be- 
ing positive and all other pseudopotentials being zero, 
and again the Hilbert space is restricted to spin-polarized 
sectors (note that for L = 1, S = 3/2 and M = 3, we 
have a dimensional vector space of wave functions, so 
no corresponding L = 1, S = 3/2 and M = 3 pseu- 
dopotcntial can occur). In these examples we observe 
that the pscudopotcntials are in fact diagonal in the L, S 
sectors. This is a direct consequence of the interaction 
being, respectively, rotationally and spin-rotationally in- 
variant. This feature will remain present in our applica- 
tion of pseudopotentials to the SSG wave function, and 
so at this point we shall drop the repeated indices and 
denote diagonal pseudopotentials by V^g = Vl [ s . LS . 

The bosonic NASS state, like the Moore-Read state, 
corresponds to a spin-polarized three-body contact in- 
teraction, however the Hilbert space now includes addi- 
tional spin sectors (not just the spin-polarized sector). 
In the language of pseudopotentials, the bosonic NASS 
state is the highest density zero energy ground state of a 
Hamiltonian with the positive, diagonal pseudopotential 
V^ 3 ^ . Motivated by the generalization of the Moore- 
Read Hamiltonian to the gaffnian Hamiltonian, our pro- 



posal for the SSG Hamiltonian is to keep V^ 3 3 ^ 2 , V^ 3 , 2 
and VJj 3 -^ positive. The proposed Hamiltonian for the 
SSG wave function is thus expressed as 

tf SSG = |0,3/2)V/ 3 3/2 (0,3/2| 
+ |2,3/2)V/ 2 3 3/2 (2,3/2| 
+ |0,1/2)^ 3 1/2 (0,1/2|. (3) 

From this pseudopotential Hamiltonian, we can al- 
ready infer some properties of its zero energy eigenstates. 
Because the three-body interaction in the 5 = 3/2 chan- 
nel is identical to the Hamiltonian generating the (polar- 
ized) gaffnian, the zero energy eigenstates will vanish as 
at least a third power when three particles with the same 
spin are brought to the same point. In addition, when 
three particles have overall spin S = 1/2, the wave func- 
tions vanish at least quadratically. These are indeed the 
vanishing properties consistent with the CFT description 
of the SSG as we shall describe in Sec. IIVI 

It is worth stating here that the SSG wave function 
likely also corresponds to the ground state of other, more 
complicated, local Hamiltonians involving M-body terms 
with M > 3, following the line of reasoning discussed in 
Ref . ill. 



B. Spin-Dressed Squeezing Algorithm 

There is another way to characterize the 'spin-singlet 
gaffnian', apart from by means of the Hamiltonian we 
introduced above. This method is the so-called 'squeez- 
ing method'. The idea underlying the method is that 
for many model states, the wave functions have a large 
number of zero coefficients if expressed in terms of the 
space of all possible monomials. These zeros are closely 
related to the vanishing properties of the wave functions. 
The method was pioneered in Ref. [29|. The spinful case, 
which we will employ here, was described in great detail 
in Ref. H 

We consider quantum Hall states on the sphere^, in 
the presence of iV$ flux quanta. In the lowest Lan- 
dau level, this construction gives rise to iV$ + 1 or- 
bitals, whose occupation numbers will be denoted by 
(no,ni, . . . jUn^). These orbitals have angular momen- 
tum (A$/2,iV$/2 - 1, ...,-N<s,/2). We will use the 
orbital occupation numbers to label the states in the 
Hilbcrt space. Because total angular momentum is a 
good quantum number, we can split the Hilbert space 
into sectors with different values of total L z , and we will 
always consider the number of particles N to be fixed. 
For spinless fermions, the Pauli-principle specifies that 
the occupation of any orbital ni can be zero or one, while 
for bosons, there is no constraint. 

To obtain all states in a particular L z sector, we divide 
the N particles over the orbitals, such that one obtains 
the correct value of total L z , and such that the particles 
are 'desqucczed' as much as possible. This means that 
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for an even number of bosons in the L z = sector, one 
considers the occupation (N/2, 0, . . . , 0, N/2). To obtain 
the other states in this sector, one generates all possible 
pairwise rearrangements of particles occupying orbitals, 
keeping the total L z fixed. In the general case, this en- 
tails to (. . • , fii, . . . , tij—i, rij, . . .) transforming into 
(. . . , rii — 1, rii + i + 1, . . . , rij-i + 1, rij — 1, . . .), where all 
the other occupation numbers remain unchanged. In this 
way, one obtains the full Hilbert space in each total L z 
sector. 

As we alluded to above, the wave functions of many 
model states have a large number of zero components, 
if expressed in the Hilbert space described above. For 
instance, if one considers the v — 1/2 bosonic Laugh- 
lin state, it suffices to construct a reduced Hilbert space, 
by starting to squeeze from the following, so-called root 
configuration, (1, 0, 1, 0, . . . , 0, 1, 0, 1), instead of the com- 
pletely desqueezed configuration (A/2,0, ... ,0, N/2). 
All states in the full Hilbert space which do not appear 
in the reduced Hilbert space have zero coefficient in the 
Laughlin state. Moreover, it turns out that there is only 
one L = state one can construct in the reduced Hilbert 
space. Thus, to obtain the Laughlin state for a certain 
number of particles, one constructs the reduced Hilbert 
space from the root configuration (1, 0, 1, 0, . . . , 0, 1, 0, 1), 
and demands that L + on a general state in this reduced 
Hilbert space gives zero. This procedure gives a set of 
equations for the coefficients, whose solution gives the 
Laughlin state, expressed in the monomial basis on the 
sphere. In a similar way, many other model states can be 
obtained. For instance, to obtain the level-fc Read-Rezayi 
states^, one only has to change the root configuration to 
(k, 0, k, . . . , k, 0, k). Interestingly, the root configurations 
correspond to the 'thin-torus' limit of the states^ - — 

To describe spinful wave functions, the squeezing 
method was augmented in Ref. l30l . The idea is to 
start with a root configuration, for which one at first 
ignores the spin degrees of freedom. To construct the 
reduced Hilbert space, one creates a set of orbital oc- 
cupations, by squeezing in all possible ways. Finally, 
one assigns spin to all the particles, in all possible ways. 
We will be interested in spin-singlet states, for which to- 
tal S is a good quantum number. So to construct the 
ground state, we work in the total S z = sector. Af- 
ter distributing the spin in all possible ways for each 
configuration, we have constructed the reduced Hilbert 
space. The spinful orbital occupations are now denoted 
by (»o,t) n <U, ni,t) "1,4.! • ■ ■ ' n JV #1 f, «iv*,4.) . 

We still have to specify how to obtain the correct wave 
function in the reduced Hilbert space. In general, we 
will be interested in ground state wave functions which 
have L = S = 0, so to obtain those, we work in the 
L z = S z = subscctors, and demand that the action of 
L + and S + on the state gives zero. It turns out that this 
is in general not enough to completely specify the wave 
functions. In addition, one has to demand that certain 
states in the reduced Hilbert space have zero coefficient. 
The SSG is an example where this happens, as we shall 



presently describe. 

Let us now discuss how the squeezing algorithm ap- 
plies to the construction of the SSG Hilbert space. We 
have already seen, in the construction of the Hamilto- 
nian in Eq. [3J that the zero energy eigenstates of the 
SSG wave function vanish as at least a third power 
when three particles with the same spin are brought 
to the same point and, in addition, when three parti- 
cles have overall spin S = 1/2, the wave functions van- 
ish at least quadratically. These observations motivate 
the use of the root configurations based on the pat- 
tern (2, 0, 2, 0, 0, 2, 0, 2, 0, 0, ... , 0, 0, 2, 0, 2), because like 
the Moore-Read state, the resulting wave function will 
vanish at least quadratically when three particles are at 
the same location. These root configurations have the 
property that two neighbouring orbitals can maximally 
be occupied by two particles, while five consecutive or- 
bitals can be occupied by at most four particles. We 
note in advance that these root configurations satisfy the 
spinful generalized Pauli principle, which we discuss in 
the next subsection. 

To construct the ground state wave function, we start 
with the root configuration described above, and then 
construct the reduced Hilbert space. To obtain the 
ground state, we then demand that both L + and S + 
act to give zero on all states in the reduced Hilbert. This 
construction, however, only guarantees that the obtained 
wave function vanishes as a second power when three 
particles of the same spin are coincident; but we really 
wanted the wave function to vanish as a third power in 
that case. We therefore enforce an additional constraint: 
That all basis states in the reduced Hilbert space have 
zero coefficient, on the condition that either one has for 
spin-up particles, no,f = 2 and ri2,t > 0, or similarly, 
for spin-down particles, one has n.04 = 2 and 7124 > 0. 
With that additional constraint, the wave function now 
vanishes as a third power when three particles are coin- 
cident with total S = 3/2. 



C. Spinful Generalized Pauli Principle 

A further method to characterize the zero energy space 
of the SSG Hamiltonian is to use the generalized Pauli 
principle for quantum Hall wave functions^ The gener- 
alized Pauli principle is most readily expressed in terms 
of partitions (or occupied orbitals). A partition A is de- 
fined to be an ordered set of N integers, {Ai, . . . , Ajv}, 
where Ai > Aj+i. 

To describe the Hilbert space of spinless particles 
in flux N$>, we restrict the integers to the set Ai £ 
{0,1,..., iV$}. We then have the following relation be- 
tween the orbital occupation numbers no, rii, . . . , Hjv s , 
which were introduced in the previous subsection, and 
the Aj forming the partition A: Namely, n^r s is the num- 
ber of i such that Ai = N^ etc. Thus in general, rij is the 
number of i such that Ai = j. 

To characterize the quasiholes of the (k, r) clustered 
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states^ 9 - we introduce the notion of (k, r)-admissiblc par- 
titions, which are partitions obeying the following condi- 
tion for all i: 

K - K+k > r. (4) 

For given fixed values of N and iV$, one identifies 
the number of admissible partitions with the number of 
quasihole states in the corresponding spectrum. For in- 
stance, the r = 2 series corresponds to the level— k Read- 
Rezayi states, the Laughlin state (k = 1) and the Moore- 
Read state (k = 2) being the simplest cases. The gaffnian 
wave function is associated with the (fc = 2, r = 3) gen- 
eralized Pauli principle. 

In order to generalize this method to the spinful case, 
we follow the argument presented in Ref. [2ft First, the 
partition is replaced by a spinful partition that mixes mo- 
mentum and spin. A spinful partition (A, a) is specified 
by N integers, {Ai, . . . , Ajy}, and now in addition, a set 
of spin indices {a±, . . . &n}, where <jj G { — 1,1}. If ei- 
ther the condition A; > Ai + i or A^ = A; + i and Oi > ai + \ 
holds, we say that (A,cr) constitutes a spinful partition. 

The spinful generalization of the (k, r) admissible par- 
titions is given by the following conditions, which have 
to hold for all i: 

A 4 - Xi+k > r, 
or A 4 - Xi+k = r - 1 an d <Ji < <Ji+k- (5) 

In Ref. [25| it is shown that the densest bosonic state 
that leads to a spinful (fc, r) admissible possible corre- 
sponds to a filling factor v = 2/c/(2r — 1) and a shift 
S = r. Note that on the sphere, the relation between the 
filling factor v and shift 6 is 

N 

It has been proposed^ 5 , that non-symmetric Jack poly- 
nomials can represent spin-singlet states that are the spin 
generalization of the spinless clustered states. In that 
case, the root partition is chosen to be a spinful (fc,r) 
admissible partition. The number of quasihole states can 
be obtained in a similar way to the spinless case, count- 
ing the number of admissible partitions. Indeed, these 
properties have been checked for the Halperin spin sin- 
glet states (k = 1, r) and the NASS state (k = 2, r = 2). 

In our application of the spinful generalized Pauli prin- 
ciple to the SSG, we employ the (k = 2, r = 3) Pauli prin- 
ciple (as for the spin-polarized gaffnian), but now using 
the spinful partitions. 

III. NUMERICAL RESULTS 

In this section we shall present numerical evidence 
demonstrating that the zero-energy states of the local 
Hamiltonian presented in Sec. HH Eq. [3l are correctly 
reproduced by both the squeezing algorithm and by the 



generalized Pauli principle which were also described in 
Sec. [II] We shall present evidence derived from both the 
quasihole spectrum and particle entanglement spectrum. 

A. Quasihole Spectrum 

The quasihole spectrum is determined by numerical 
exact diagonalization of the Hamiltonian for finite-sized 
systems in the sphere geometry. This is done for a variety 
of system sizes, N, and for a variety of fluxes 7V$. Eigcn- 
states are labelled by the quantum numbers L z and S z 
and fall into (L, 5)-multiplets. (Note that these quantum 
numbers are completely separate from the L and S de- 
scribing the pseudopotentials unless one has only M par- 
ticles with an M-body interaction.) Given constraints on 
the dimension of the spinful Hilbert space, we were able 
to study systems of up to N = 12 and N$ = 12. An 
example of such a spectrum is plotted in Fig. [TJ 

For comparison, we generated the same Hilbert space 
via the spin-dressed squeezing and generalized Pauli 
algorithms. For the spin-dressed squeezing algorithm 
we started from a root partition (2,0,2,0,0), e.g. 
(2,0,2,0,0,2,0,2,0,0,2,0,2) for N = 12, and we ap- 
plied the procedure described in Sec. lIIBI to generate the 
Hilbert space. The ground state wave functions, corre- 
sponding to the densest root configurations with L z = 0, 
occur at filling v = 4/5, and shift 6 = 3. For the gener- 
alized Pauli algorithm, we specified a spinful admissible 
partition (k = 2, r = 3). 

Our key observations are as follows: First, in the spec- 
trum of the Hamiltonian there is a unique zero energy 
state occurring only in the L = 0, S = sector for N = 4 
at JV* = 2, for N = 8 at N<s, = 7 (see Fig. [TJ and for 
N = 12 at Ng, = 12, which all correspond to v = 4/5 and 
<5 = 3, consistent with both the squeezing and general- 
ized Pauli approaches; second, for the quasihole spectrum 
generated for N between 3 and 12 and for up to 20, 
we have checked that the counting of zero energy states 
in each L z and S z sector precisely matches the counting 
predicted by both the spin-dressed squeezing algorithm 
and the generalized Pauli principle. 

We have further determined for N = 4 and N = 8 that 
the ground state monomial expansion of the Hilbert space 
generated by exact diagonalization of Eq. [3] precisely 
matches the ground state monomial expansion generated 
by the squeezing algorithm as described in Sec. Ill Bl In 
addition we have found that the ground state monomial 
expansion also precisely matches that of the proposed an- 
alytic form of the ground state SSG wave function (we 
shall discuss the analytical form of the ground state wave 
function in Sec. IVl). 

B. Particle Entanglement Spectrum 

The so-called particle entanglement spectrum (PES) is 
determined from the reduced density matrix of a subsys- 
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FIG. 1: The 5 = sector of the full spectrum of the 
Hamiltonian -£fc;SG> denned in Eq. [21 for N = 8 and 
A$ = 7, obtained by numerical exact diagonalization in 
the sphere geometry. The inset zooms in on the bottom 
left corner of the spectrum and shows the unique zero 
energy state located in the L = 0, 5 = sector. 



ical calculations here, is what is the size of the gap for 
quasi-particle excitations in the thermodynamic limit? 
For finite-sized systems, looking at Fig. [TJ the gap is 
clearly finite (as it is in the gaffnian case) . There is good 
reason, however, to argue that the gap will in fact van- 
ish in the thermodynamic limit, and hence the SSG will 
be gapless and the wave function therefore compressible 
(see Ref. [40h . That argument is based on the fact that 
the CFT corresponding to the SSG is non-unitary, as we 
shall presently discuss. 



IV. CONFORMAL FIELD THEORY 

Let us now discuss the guiding principles of the deriva- 
tion of the CFT for the SSG wave function (see also re- 
lated considerations in Ref. l25l ). We shall expand further 
on the discussion in appendix E] For an introduction on 
CFT, we refer to the book, Ref.|4l|, and the seminal paper 
by Belavin, Polyakov and Zamolodchikov (BPZ), Refl42i. 



tem that results from the partition of the whole system 
into two (or more) parts A and Bi 31 i 32 During this par- 
titioning we also keep the overall geometry unchanged. 
The reduced density matrix pa is given in terms of the 
full density matrix p = \ip) (ip\ by tracing out the Nb 
particles in the B partition: pa = Ttb(p). The PES 
arises from diagonalizing pa and then classifying the re- 
sulting eigenstates according to the symmetries of the 
problem, in our case by L z and S z . 

For model states for topological phases such as the 
Laughlin and Moore-Read states, it has been observed 
that there is a characteristic PES: The number of non- 
zero eigenvalues for pa is identical to the number of quasi- 
hole states for a system with identical geometry but only 
Na particles i 31 ' 32 The number is usually exponentially 
lower than the dimension of pa, or equivalently, there is 
an infinite so-called 'entanglement gap' to the remaining 
eigenvalues. 

For the SSG wave function, our key observations are 
as follows: We find evidence of an infinite entanglement 
gap for N = 4 at A$ = 2, for N = 8 at iV$ = 7 and for 
N = f 2 at = 12, which we associate with our model 
ground state wave function for the SSG; and we find that 
the counting predicted for the quasihole excitations up to 
N = 12 precisely matches the counting predicted by the 
generalized Pauli principle in each 5 and L sector. 

With the evidence presented in this section in mind, 
we conclude that the proposed SSG Hamiltonian in Eq.[3] 
generates the correct zero energy Hilbcrt space for the 
SSG wave function, at least for system sizes up to N = 
12. Based on this evidence, we also expect our conclusion 
to hold for any other system size. Thus we have a Hamil- 
tonian which describes the ground state and excitation 
spectrum of the Hilbcrt space. 

An important question, left unanswered by our numcr- 



A. Coset constructions for minimal models 

The minimal models introduced by Belavin, Polyakov 
and Zamolodchikov can be written in terms of coset mod- 
els of the su(2)k WZW models In the unitary case, one 
has 



M(k + l,k + 2) 



SU{2) 1 X SM(2) fc _! 

su{2) k 



(6) 



For k = 2, this gives the Ising CFT, M(3,4). For inte- 
ger k, the coset theories are unitary, a property which is 
inherited from the su(2)k WZW model, which is unitary 
for k integer. 

General minimal models are labeled M(p',p) for arbi- 
trary (non-negative) co-prime integers (p',p). The mini- 
mal model is unitary only if \p — p' \ = 1 , otherwise it is 
non-unitary. Nevertheless, there exists a coset descrip- 
tion of non-unitary minimal models in terms of WZW 
models at fractional leve h 44 ' 45 In particular, the param- 
eter k is given by k = 3 p~p, p ■ The non- unitary minimal 
model A4(3,5), featuring in the CFT description of the 
gaffnian wave function, corresponds to k = —1/2, while 
the Yang-Lee model M{2, 5) has k = —4/3. 

The central charge of the minimal models is given by 
c (p\p) = 1 — 6 ^ P pp /P ' > ' The primary fields of the mini- 
mal models are labeled by integers (r, s), which take the 
values 1 < r < p' and 1 < s < p. The labels (r, s) 
and (p' — r, p — s) correspond to the same primary field. 
Finally, the conformal dimensions of the fields are given 
by 



h(r,s) 



(rp — sp')" — (p' — pY 
4pp' 



(7) 
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B. Gepner parafermions in terms of minimal 
models 



The non-abelian part of the CFT describing the 
NASS state, was originally written in terms of Gepner 
parafermions 4 ^, which can be expressed in terms of the 
coset 



GPf(«u(3) fc ) = 



su(3) k 



u(l) 2fc u(l) 6 fc 



For our present purposes, we shall focus on k = 2, in 
which case the above coset is equivalcntly written as 



GPf(sw(3) 2 



su(2)i x su(2)i x su(2)i 
su(2) 3 



(8) 



(for arbitrary su(3)k parafermions, su(2) has to be re- 
placed with su(k) in the above). If one multiplies the 
numerator and denominator of the coset above by su(2)2, 
one can factorize the coset into the (semi-direct) product 
of the minimal models M(3, 4) and .M(4, 5). 



GPf(su(3) 2 ) 



sw(2)i x su(2)i x su(2)i x su(2) 2 

su(2) 2 x su(2) 3 
J\A (3,4) x A4 (4, 5). (9) 



Let us be more precise about this correspondence. The 
direct product .M(3,4) x A4(4, 5) does not correspond 
to the su(3)2 Gepner parafermions. Instead, one has to 
consider a so-called different modular invarian t 47 ' 48 , con- 
structed from the fields present in the direct product. To 
do that, we will follow the logic presented in Refs. l49T - [5ll 
This amounts to identifying a boson in the CFT, which 
is then said to be 'condensed' (or added to the chiral al- 
gebra). To identify a suitable boson, we give the fields of 
the models A^(3,4) and .M(4, 5) in table HJJFusion rules 
for the minimal models are given in RcfsTfill andl42h 



M(3,A) 


h(r,s) 


s = 1 2 3 


r = 1 


1/16 1/2 


2 


1/2 1/16 



h(r, 


s) 


s = 1 2 3 4 


r = 


= 1 


1/10 3/5 3/2 


r = 


- 2 


7/16 3/80 3/80 7/16 


r - 


= 3 


3/2 3/5 1/10 



TABLE I: Kac table of conformal weights for the 
unitary minimal models A4(3,4) and A / I(4,5). 

We will label the fields in the product theory in 
terms of the conformal dimensions (or conformal weights) 
of the contributing fields of the original CFTs. We 
find that there is indeed a bosonic field (i.e., a field 
with integer scaling dimension) in the product, namely 
(1/2,3/2). The condensation picture amounts to the 
following procedure 4 ^: Particles which can be obtained 
from one another by fusion of the boson, are 'identified' 
(in the same way as the boson itself is identified with the 
vacuum, or trivial particle); in addition, particles which 



are not mutually local with the boson, are 'confined'; 
a particle which is mutual local with the boson does not 
generate a phase factor when transported around the bo- 
son; finally, it can happen that particles 'split'. We will 
see an example of splitting below. In this section, we will 
be rather brief, more details on how one constructs the 
correct theory can be found in the appendix. 

To get started, we list the particles in the product the- 
ory which are mutually local with this boson, namely 

e= (0,0) = (1/2,3/2), 

(-^)„ = (1/16,3/80), 



(0,1/10) = (1/2, 3/5), 



(h a = (1/16, 7/16), 



3 
5 



(0,3/2) = (1/2,0), 
(0,3/5) = (1/2, 1/10). 



The equalities within one line signify that the particles 
corresponding to the two fields are identified. As an 
example, one obtains (1/2,3/5) from (0,1/10) by fu- 
sion with the boson, namely (0,1/10) x (1/2,3/2) = 
(1/2,3/5). We also note that the particles (|) a and 
(yq)o have to split into two particles, thereby giving the 
full Z3 symmetry present in the sw(3)2 parafermion the- 
ory. That this splitting is necessary can be seen by 
taking for instance the fusion product (i) a x (i) a = 
(0,0) + (1/2,0) + (0,3/2) + (1/2,3/2). Because both 
(0,0) and (1/2,3/2) correspond to the identity, split- 
ting of the field (\) a is necessary. A similar argument 
applies for {jfi)a, and one can convince oneself that in 
the end, one indeed obtains the correct fusion rules of 
the sm(3)2 parafermion theory. We conclude this sec- 
tion by mentioning that the example of the equivalence 
GPf(sw(3) 2 ) « M(3,4) k M(4,5) featured prominently 
in a paper relating the Moore-Read and NASS states^ 



C. CFT for the spin-singlet gaffnian 

The spin-singlet gaffnian that we aim to construct 
should have the same vanishing properties as the gaffnian 
when considered as a function of only one spin species of 
particles. Motivated by the product of minimal models 
describing the su(3)2 Gepner parafermions, we take as a 
starting point the coset 



su(2)i x su(2)\ x su(2)fc 
sw(2) fe+2 

with k some fraction. Multiplying numerator and de- 
nominator by su(2)k+i, we can write 

su(2)i x su(2)k su(2)i x su(2)k+i 
— X 



su{2) 



k+l 



su{2) 



fc+2 
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The CFT describing the gaffnian is the non-unitary min- 
imal model -M(3, 5). To retain 5) in the above con- 
struction, we will choose k = —1/2. With this choice, our 
coset takes the form 5) x VW(5, 7). In analogy with 

the Gepner parafermion case above, we look for a field 
with integer scaling dimension in the product theory, and 
condense it. The field content of the models is given in 
table im (F usion rules for the minimal models are given 
in Rcfs. S3 and |H) 



M(3,5) 


h(r, a) 


s = 1 2 


3 4 


r = 1 
2 


-1/20 
3/4 1/5 


1/5 3/4 
-1/20 



M(5,7) 



h(r, 


s) 


s = 1 


2 


3 


4 


5 


6 


r — 


1 





1/28 


3/7 


33/28 


16/7 


15/4 


r = 


2 


11/20 


3/35 


-3/140 


8/35 


117/140 


9/5 


r — 


3 


9/5 


117/140 


8/35 


-3/140 


3/35 


11/20 


r = 


4 


15/4 


16/7 


33/28 


3/7 


1/28 






TABLE II: Kac table of conformal weights for the 
unitary minimal models A4(3, 5) and .M(5, 7). 



It turns out that the only boson that we can condense 
(or add to the chiral algebra), is the field (9/5,1/5), 
which curiously has non-trivial fusion rules. The only 
fields that survive this condensation, are found to corre- 
spond to product of fields of the first column of the Kac 
table for A^(3,5) and the fields of the first row of the 
Kac table for Ai(5,7). In addition, none of these prod- 
uct fields are split. It follows that the resulting theory 
can be thought of as the product of a non-unitary semion 
theory (two fields with Zi fusion rules and dimensions 
and 3/4) and a theory of six fields, satisfying su(2)^ fu- 
sion rules. This second theory is also non-unitary, and 
the scaling dimensions are again simply read off form the 
Kac table, (0,1/28,3/7,33/28,16/7,15/4). We note that 
this theory is modular (inherited via the coset construc- 
tion). Here, we were rather brief in our description of 
the construction of the theory M(3, 5) X M(5, 7). In ap- 
pendix IA 21 we will give the details of the construction. 

Before we start with the construction of the electron 
and quasihole operators, in table IIIII we give the Kac 
table of the fields present in the CFT we constructed, 
Ai(3, 5) k A4(5, 7), which will constitute the non-abelian 
part of the CFT describing the SSG. The full CFT 
also contains the u(l) vertex operators (see for instance 
Ref. I53T ) . It is important to note that one has to be care- 
ful in determining the scaling dimensions of the fields 
in the product theory. The dimensions obtained by 
simply adding the scaling dimensions of the constituent 
fields can in fact correspond to the scaling dimensions 
of descendant fields, differing from the scaling dimen- 



■M(3,5) k M(5,7) 





3 = 1 2 


3 4 5 


i = 
1 


1/28 3/7 
3/4 11/14 5/28 


5/28 2/7 3/4 
-1/14 1/28 1/2 



TABLE III: Kac table of conformal weights of the CFT 
describing the non-abelian part of the spin gaffnian. 



sions of the primaries by integers. In particular, the 
field (15/4,0) is identified with (9/5,1/5) x (15/4,0) = 
(11/20, 1/5), which has scaling dimension 3/4. Similarly, 
(9/5, 1/5) x (15/4, 3/4) = (11/20, -1/20), which has scal- 
ing dimension 1/2. 

The fields occupying the corners of the Kac table, 
namely (0, 0), (0, 3/4), (15/4, 0), (15/4, 3/4), are special, 
because they are simple currents. A simple current is a 
field, which when fused with any other field, always gives 
a single field as the result. Therefore, the particles cor- 
responding to simple currents are abelian. To construct 
the 'electron' operator, one is therefore only allowed to 
use one of these four fields. Before we get started, we 
note that the abelian sector of the theory at hand has 
less symmetry than the one describing the NASS state, 
which has three primary fields all with the same confor- 
mal weight. 

We will now work with the assumption that we do have 
a Z2 symmetry between the fields (15/4, 0) and (0, 3/4), 
(i.e., the fields with labels = (0,5), (1,0)) which 

we need, if the correlator is to describe a spin-singlet 
state (in the appendix, we will give evidence supporting 
this statement). In addition, we want the wave function 
to look like a gaffnian when viewed as a wave function 
for either spin-up or spin-down particles alone. We will 
introduce the following notation for the fields 

1 = (0,0), Vt = (0,3/4), 

V^ = (15/4,0), i> n = (15/4,3/4). (10) 

The fusion rules of these fields read ipf x ip^ = ip^ x ijj^ = 

■0ti x i/>n = 1 and x i>x = V>n- 

We will continue with adding the appropriate vertex 
operators, giving charge and spin to the constituent par- 
ticles, which is the usual procedure (see for instance 
Ref. 1531 ) . The guiding principle will be to construct elec- 
tron operators, which give rise to the vanishing properties 
we want. The ansatz for the operators is 

V i (z i ,)=1>M)e ia *'- m '(zi), (11) 

where a and /? are constants to be determined. The 
fields <p c and (j> s are u(l) compactified bosons. As we 
already pointed out, when we bring several up-particles 
together, we demand that the wave function behaves in 
the same way as the gaffnian wave function. The 'elec- 
tron' operator in the CFT description of the gaffnian 
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reads V{z) = ip{z)e % ^i^{z), with t/)(z) a field with con- 
formal dimension 3/4. Equivalences with this polarized 
case gives us the constraint a 2 + (3 2 — 3/2. 



Moreover, brining an spin-up and spin-down particle 
together should not lead to a divergence. Using the fusion 
rule Vf x V\. = '0ti e2ia ^ c gives rise to a factor (z-j- — 
z^+^-P in the operator product expansion. To avoid 
a pole, one obtains the constraint a 2 — f3 2 = p, with p 
an integer greater or equal to one. Picking the minimal 
choice p = 1, one finds that a = \/5/4 and (3 = 1/2. This 
value for j3 is in fact equal to the value this parameter 
takes in case of the NASS state. 

With the electron operators V t = ^ e iV^^+l/2i<t>. 
and Vj. = ^eiyfi/Ztc-l/Wa m p i acej we conclude that 
the corresponding wave function does not vanish when 
any two particles come together. It vanishes as three 
powers when three up (or three down) particles come 
together. In addition, when two up and one down parti- 
cle come together, the wave function generically vanishes 
quadratically. We then assume that one can use logic 
similar to the NASS case, and thus, exploiting the SU (2) 
symmetry^, one finds that the wave function vanishes as 
three powers when three particles in the 5 = 3/2 channel 
are coincident, and quadratically when three particles in 
the S =1/2 channel are coincident. 



The filling fraction of the wave function is given by v = 
4/5. The scaling dimension h e — 3/4 + 3/4 = 3/2 gives 
rise to a shift S = 3, the same as for the gaffnian. We will 
defer the construction of the quasihole operators to the 
appendix, but mention here that they are constrained by 
the fact that the wave functions for the electrons should 
not have poles, even in the presence of quasiholes. For 
N particles the ground state SSG wave function can be 
written as a conformal block containing equal numbers 



of spin-up and spin-down electron fields: 

*SSG = (M4)---M*\)M*$ +1 )"-M*n)) ■ 

Finally, we would like to mention that if we would 
have assumed that the field V>t4 would have had con- 
formal dimension 3/2 instead of 1/2 (which is not consis- 
tent, because of the field identification above), one would 
have found a = = \/3/4. The resulting wave function 
would then factorize as gaffnian(up) x gaffnian(down). 
The CFT description for such a wave function should 
be M(3,S) x M(3,5), which is at odds with the CFT 
description we have used here. 

V. GROUND STATE WAVE FUNCTION 

In this section we shall present an explicit construction 
for the ground state SSG wave function for the bosonic 
case. Our construction is conjectured with the view to 
satisfying the vanishing properties as well as the con- 
straints of filling factor v = 4/5 and shift <5 = 3, arising 
from the CFT considerations put forward in Sec. IIVI In 
addition, the wave function must describe a spin singlet 
state. That condition specifies the requirement to satisfy 
the Fock cyclic symmetry conditions (see Ref. I54T ). 

Before we describe how to construct the SSG, it will 
be useful to motivate our methodology by briefly review- 
ing the construction of the spin-polarised gaffnian wave 
function^ For the gaffnian, the CFT requirement im- 
posed on the vanishing properties arc that for any three 
particles coincident the wave function must vanish as 
three powers and that the wave function must not van- 
ish for any two particles coincident (the CFT describing 
the gaffnian is the non-unitary minimal model -M(3, 5)). 
Further, it was determined that in its bosonic form, the 
gaffnian occurs at filling 2/3 and 5 = 3. 

To give the expression of the gaffnian wave function 
that can be generalized to the SSG, we divide the parti- 
cles into two groups A and B of equal size. The gaffnian 
wave function can then be written as 



^gaffnian 



= S 



n o 

i<j£A 



4f n << 

<jeB 



-4? II (^"^)Per 

i£A,j£B 



where S represents a symmctrization operation over all 
N particle coordinates, and Per[My] denotes the 'per- 
manent' of a matrix M whose elements are in this case 
given by Mjj = (zf — z^) -1 . On its own, this perma- 
nent factor contains a certain pattern of poles. When 
placed within the full construction, these poles conspire 
to ensure that the overall wave function does not vanish 



as two particles become coincident. Thus the polyno- 
mial vanishes only when three particles are coincident. 
The vanishing power can be tuned by adjusting the ex- 
ponents of each Jastrow-type factor in the construction. 
For example, W^^A^t ~ z f) 2 cou ld be adjusted to, 
say, rii^eA^i 4 — z f) A (and similar for the particles in 
group B) to alter the vanishing power from 3 to 5 in this 
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example. An entirely non-trivial step is to determine 
whether or not these vanishing properties are retained 
once the overall symmetrization operation has been com- 
pleted. For the gaffnian, it was found that the vanishing 
properties are retained. 

Now let us turn to the construction of the ground state 
wave function for the SSG. We have determined that the 
following construction, given in Eq. 1121 gives a wave func- 
tion at filling factor v = 4/5 and 5 = 3 and it satisfies the 
CFT vanishing constraints for N = 4, 8, 12. We have rea- 



son to believe that it will work for all other values of N 
(i.e., 16, 20, etc.), although an explicit check is not pos- 
sible. Thus, we conclude that it must be proportional to 
the conformal block giving the SSG ground state. In the 
construction, we have used a clustering principle: The 
particle coordinates are first divided into equal sets of 
spin up and spin down, and then each of these sets is 
further subdivided into two equal subsets labelled by ei- 
ther A or B, giving four sets in total (A A I, B 'f, 
B I). With the label A or B on its own, we refer to the 
two sub-sets with either spin direction. 



SSG 




n o 

i<j€af 



at 



n 



ai\2 



at 



n w 



J ) Per 



j 

(12)" 



where, as above, Per[Afij] denotes the 'permanent' of 
a matrix M whose elements are in this case given by 
My = (zf — zf) -1 and where Yg—o is the Young operator 
for a spin singlet representation of the symmetric group — 
this operation is required in order to guarantee that the 
wave function satisfies the correct Fock cyclic symmetry 
conditions for a spin singlet stated 

Ys-o = S t t S i i A t i ...A t i , 

1" N/2 1 '"*N/2 11 N/2 N/2 

that is, the Young operator for a spin singlet represen- 
tation corresponds to the operation of antisymmctrizing 
over ordered pairs of spin-up and spin-down coordinates, 
followed by symmetrizating over all spin-down and then 
all spin-up coordinates. Crucially, it is important to note 
that the wave function does not completely vanish when 
Ys=o is applied! 

In Sec. IIIII we provided strong evidence showing that 
that the spectrum of the Hamiltonian for the SSG con- 
tains a unique zero energy state with L = and S = 0, 
for precisely the values of particle number and flux cor- 
responding to filling v = 4/5 and 5 = 3. Therefore we 
conclude that this trial ground state wave function is 
unique, and further, it is precisely the highest density 
ground state of the SSG Hamiltonian proposed in Eq. [3] 

VI. DISCUSSION 

To summarize our findings, we have presented evidence 
to demonstrate that the proposed SSG wave function sat- 
isfies many of the 'ingredients of an ideal theory of a 
FQHE state': A local Hamiltonian, a relatively simple 
analytic expression for the wave function, at least for the 
ground state, and a correspondence of that wave func- 
tion to a 2D rational CFT. We have verified the quasi- 
hole spectrum of the SSG Hamiltonian, checking against 
both the spin dressed squeezing algorithm and the spin- 
ful version of the generalized Pauli principle. Indeed, 



this study is an interesting test case for the application 
of such techniques. 

Concerning quasiholcs, although we have made 
progress on the CFT description of the quasihole op- 
erators (see appendix [A| , we have not yet determined 
simple, analytic expressions for the corresponding quasi- 
hole wave functions. For the spin-polarized gaffnian, the 
quasihole wave functions have been constructed^; the 
main stumbling block here is the additional complexity 
due to the spin degree of freedom. 

In Sec. [H] we described how to generate the SSG 
wave functions by means of a local Hamiltonian writ- 
ten in terms of spin-dependent pseudopotentials. An 
interesting question left unaddressed is whether other 
states could be constructed with faster vanishing prop- 
erties than the SSG. For instance, in the spin-polarised 
case, adding the next highest L pseudopotential to the 
gaffnian Hamiltonian is known to produce the Haffnian 
Hamiltonian^ Might we be able to generate a 'spin- 
Haffnian' state with a Hamiltonian containing positive 
^0 3/2' ^23/2 ' ^01/2 anc ^ now > m addition, V^ 3 ^ 2 (and 
possibly V^ 1 i^)7 Wave functions of this type were also 
considered in Ref. HjJ Presently we lack a correspond- 
ing CFT description with which to conduct the same 
checks as for the SSG. It is also unclear at the present 
time whether that Hilbert space could identically be con- 
structed by a squeezing algorithm or generalized Pauli 
principle approach. 

Concerning CFT coset constructions, it is notewor- 
thy that we were able to derive a self-consistent CFT 
from the product of two non-unitary minimal models, by 
constructing a non-diagonal modular invariant. Specif- 
ically, the non-unitary models in the product can be 
thought of as cosets of fractional level affinc Lie algebras. 
The resulting CFT corresponds to the non-unitary coset 
su{2)\ x su(2)i x su(2)k~2/su(2)k, with k fractional. In 
the appendix IA 21 we deal with a more general case. 

In general, for any unitary coset construction with sim- 
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pie current extensions there is a proven procedure to gen- 
erate the CFT i 50 i 51 For non-unitary cosets, no such gen- 
eral procedure exists, and each CFT must be constructed 
on a case-by-case basis4^ In the CFT construction pre- 
sented in this work, we have provided a further example 
of a case where the construction of non-diagonal modular 
invariants for non-unitary CFTs is possible. 

Another powerful tool for the analysis of the quan- 
tum Hall wave functions has been the concept of Jack 
polynomials , 25 1 29 The Jack polynomials provide a conve- 
nient basis in which to describe polynomial wave func- 
tions with precisely defined vanishing properties. Based 
on this new insight, it has been determined that states 
such as that of Moore and Read — where there exists a 
corresponding local Hamiltonian, a relatively simple an- 
alytic form for the wave function, and a CFT description 
of the state — are in fact simply special cases within a 
much broader classification of FQHE states in terms of 
these Jack polynomials. It has been further realized that 
such states are rather atypical: Out of all the Jacks, there 
exist only a handful of special cases for which all three of 
these ideal ingredients of a theory of a FQHE state are 
present. In this paper we have presented a further ex- 
ample, albeit that the CFT corresponding to the SSG is 
non-unitary implying a gapless, compressible wave func- 
tion. 
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Appendix A: More Details about the CFT 
1. Quasihole Operators 

In this appendix, we will construct the operators cor- 
responding to the quasiholcs with the smallest quantum 
numbers. From them, wc construct the set of all quasi- 
hole operators, whose number gives the torus degeneracy. 

Wc start with the field in the CFT M(3, 5) x M(5, 7), 
and add the u(l) factors (as we did in the construc- 
tion of the electron operators), in such a way that the 
quantum numbers are minimized. The fields in the 
M{3,5) k M(5, 7) theory are denoted by faj, with 
i = 0,1 and j = 0,1,2,3,4,5. The label i adds mod- 
ulo 2 under fusion, while the label j obeys su(2)s fusion 



rules. For general k, the su(2)k fusion rules read 



m-\i\(j 1 +j 2 ,k-j 1 -j 2 ) 

h x h = 3 3 ■ 

h=\3i—h\ 

The two fields with the smallest scaling dimension are 
0o,i an d 0i,4, so we start by writing the ansatz for the 
smallest quasihole operator as V q h(w) = 0o,ie m ^ c+4&< ^ s . 
Fusing this quasihole with an electron should give a wave 
function which is still analytic in the electron coordinates. 
This leads to the constraints o-y/5/4 + 6/2 — 1/2 = p 
and a-v/5/4 — b/2 = q, where p and q arc non- negative 
integers. Upon picking the minimal choice p = q = 0, 
one finds a = l/(2y/5) and b = 1/2, which corresponds 
to a quasihole with charge 1/5 and spin s z = 1/2. The 

operator reads V q h^ = 0o,ie I 2V^ c+i / 2( ^ . 

Constructing the other minimal quasihole, based on 

01,4, one obtains Vqh,i = 0i,4e l 275^ c_i ^ s , with charge 
1/5 and s z = —1/2. Starting from these 'fundamental' 
quasiholcs, one obtains all the quasihole species by suc- 
cessive fusion. 

When wc construct the other quasiholcs, wc will do 
this 'modulo the electron operators', because two parti- 
cles which can be obtained from each other by fusion of 
an 'electron' correspond to the same species of particle. 
Thus, we will first determine the operators which one can 
obtain by fusing (possibly different species of) electrons. 

Wc first introduce some notation for the labels of the 
different fields (or species of particles). The fields will 
be denoted by their three 'quantum numbers', namely 
the field of the non-abelian CFT, the charge and the 
spin. As an example, the 'electron' operators are given 
by (0o,5, 1) ^) and (01,0; 1,— h)- The 'composite' of these 
two operators is (01,5,2,0). It will be useful to con- 
sider the 'dual' operators, which read (0o,5, — 1, — |), 
(01,0,-1,1) and (01,5,-2,0). From these operators, 
we can construct particles without charge, for instance 
(00,5,1, |) X (01,0,-1, |) = (0i,5,O,l) and (0i, o ,l,-|)x 
(0o,5, — 1, —h) — (01,5,0,-1). Since fusing a quasihole 
with any of the operators above does not give us a new 
species of quasihole, it follows that the set of independent 
quasiholes can be labeled such that the charge and spin 
take the values < q < 1 and s z = 0, 1/2, respectively. 

The quasihole operators we constructed above can 
be written as (0o,i, g, |) and (0i,4, g,— \). These two 
fields are in fact to be identified, because (0i,4, \, — |) X 
(0i,5, 0, 1) = (0o,i, 5, \)- Thus, to construct all the dif- 
ferent topological sectors, it suffices to repeatedly fuse 
the quasihole (0o,i,^,^), and record the different sec- 
tors, modulo the 'electron' operators. 

This procedure leads to the following topological sec- 
tors: 
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(00,4,0,0) 

(00,2,0,0) 

(00,0,0,0) 



(01,1,1,0) 

*0,3, §,§) (01,2 
(01,3,1,0) 



3 

> 5' 2^ 



V4,|,0) 

4o,3, |,0) 



i) (01,4, f, 2) 

(01,5, |,0) (00,0, |,0) 



5' 2' 



The fields which would appear in the next column, with 
charge 1, can be identified to the fields in the first column, 
with charge 0. We thus find that the number of different 
topological sectors is 15, which is the degeneracy on the 
torus. 

It is interesting to note that these sectors can also be 
understood from the so-called 'thin-torus limit— ~—, the 
generalized Pauli principle^ 5 ., or the squeezing patterns 
we discussed in section Hi B[ 

First, the pattern which we used to define the spin- 
singlet gaffnian via squeezing, is (2,0,2,0,0), i.e. two 
neighbouring orbitals can host maximally two parti- 
cles, and three through five consecutive orbitals can 
be occupied by maximally four particles. This gives 
rise to the sectors (2,0,2,0,0,2,0,2,0,0) (five times 
on a torus), (2,0,1,1,0,2,0,1,1,0) (five times on a 
torus) and (1,1,1,1,0,1,1,1,1,0) (five times on a 
torus). Each of the fields (0o.o, 0, 0), etc. corre- 
sponds to one pattern. In the first column, we find 
(00,0,0,0) ee (2,0,2,0,0), (0o,2,O,O) ee (1,0,1,1,1) and 
(00,4,0,0) ee (0,2,0,1,1). The other fields can be 
obtained by fusing with (00,1,5,2), which brings one 
from one column to the adjacent one on the right. 
On the level of the patterns, one has to to 'hop' 
one particle one step to the right, assuming periodic 
boundary conditions. Thus, one finds (2,0,2,0,0) — s- 
(2,0,1,1,0), while in the other two cases, there are 
two options, (1,0,1,1,1) -> (2, 0,1, 1,0); (0,1, 1,1,1) 
and (0,2,0,1,1) -> (0, 1, 1, 1, 1); (0, 2, 0, 0, 2). Here, 
we have the identification (00,1,5,2) = (2,0,1,1,0), 
(00,3,1, |) = (0,1,1,1,1) and (0o,5 J5 ;, D = (0,2,0, 0,2). 
The remaining columns follow in the same way. 



2. Field content and the modular invariant of the 
product CFT 

To completely specify the conformal field theory de- 
scribing the spin-singlet gaffnian, we need to give the 
field content of the product CFT we use. In section HVl 
above, we used the condensation picture to argue which 
fields are present in the final CFT we used. We did not 
in fact go through the condensation picture in full detail: 
In particular, we glossed over the fact that the boson it- 
self has to split (because the fusion of the boson with 
itself has multiple fusion channels, including the vacuum 
and the boson itself). So, in the construction of the SSG 
CFT it is only one part of the boson which condenses, 
and going through the whole condensation procedure can 
become cumbersome. 



In this appendix, we will construct a modular invariant 
partition function of the product theory, using the boson 
condensation as a guiding principle. The first step of the 
procedure is to find a boson which can be condensed and 
to find the fields local with respect to that boson. Such 
fields will be present in the final CFT. Upon acting with 
the boson on these fields, we obtain the fields (in the 
original product theory) which are identified with those 
fields which are local with respect to the boson. In this 
way, we can construct the modular invariant partition 
functions. We refer to Refs. I47l - l49l for more details. 

We start with those CFTs corresponding to wave func- 
tions with n = 2 components, i.e. the generalizations 
of the NASS state, where the condensate particles have 
spin-1/2. The CFT describing these states consists of a 
product of two minimal models, which, for r > 2, are 
non- unitary (the case r = 2 is the NASS state). 

We already argued that for r = 3, the starting point is 
the product theory M(3, 5) x M(5, 7). For general r, we 
instead find that we need M(3, r + 2) x M(r + 2, 2r + 1) 
as the starting product theory. We note that for the 
minimal model theory M.(p',p) to be well defined, p' and 
p have to be co-prime, which means that r mod 3+^1. 
We will assume that this relation is satisfied from now 
on. 

The fields of the product theory will be labeled by 
(ri, Si; T2, S2), with the usual field identifications, com- 
ing from the field identifications present in the mini- 
mal model M(p',p). In the product theory, the field 
6 = (3,1; 1,3) always has scaling dimension h\, = 2, 
and hence is a boson. For r = 2, we find 6x6 = 
(1,1; 1,1), the identity. For r = 3, we have 6 x b — 
(1,1; 1,1) + (1,1; 1,3) + (3,1; 1,1) + (3,1; 1,3). For r > 5 
(but r mod 3 + 1 1), we find the general result 6x6 = 
(1, 1; 1, 1)+(1, 1; 1, 3) + (l, 1; 1, 5) + (3, 1; 1, l) + (3, 1; 1, 3) + 
(3,1; 1,5) + (5,1; 1,1) + (5,1; 1,3) + (5,1; 1,5). 

We now describe the field content of the CFTs we 
employ in this paper, by specifying the relevant modu- 
lar invariant partition functions. In general, these are 
given in terms of the characters of the holomorphic 
and anti-holomoriphic part of the CFT. We will denote 
these characters by Xr and \ r , respectively, where both 
r are r' are labels of the fields in the product CFT. 
The putative partition functions take the form Z(t) = 
J2 r r' ^r,r'Xr(T")x r '(T), where the M r y are non-negative 
integers denoting the multiplicities of the holomorphic- 
anti- holomorphic combination Xr-( T )Xr'(^)- The identity 
field should be present and non-degenerate, M^i = 1. 
In order that the partition function is modular invari- 
ant, the matrix M should commute with both S and T, 
the generators of modular transformations. The com- 
mutation condition can be seen to enforce that the only 
combinations XX which can occur have h — h mod 1 = 0. 

The so-called diagonal invariant, M r y = 5 r y always 
exists, but this is not the invariant we are interested in 
here. We are interested in the invariant, which contains 
4r fields, which are 'composed' out of several fields, which 
arc identified via the condensation of the boson present in 
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the original product theory. To specify the invariants, we 
have to distinguish two cases separately. We start with 
r odd. In this case, the fields in the final theory can be 
labeled by (1, s\; r2, 1), for s± = 1, 2, and r% = 1, . . . , 2r, 



for a total of 4r fields. These fields arc identified with the 
fields (q, si, r%, q), with q = 1,3, ...,r, i.e. all possible 
odd values of the labels r\ and si. In particular, the 
modular invariant partition function takes the form 



2r 2r 

|2 



Z ( T ) = E |X(l,l;r 2 ,l) + X(3,l;r 2 ,3) H ^ X(r,l;r 2 ,r) | + |X(l,2;r 2 ,l) + X(3,2;r 2 ,3) H 1" X(r,2;r 2 ,r) | • (Al) 

r 2 =l r 2 =l 

In the case that r is even, the new theory still contains the fields which are labeled by (1, 1; r%, 1), but the fields with 
labels (1, 2; r%, 1) are absent. Instead, the fields with labels (2, 2; j"2, 2) are now present. These are identified with the 
fields (q, 2; r-i, q), where q = 2,4, . . . , r is even. The partition function reads 

X(l,l;r2,l) X(3,l;7*2,3) ' ' ' ~\~ X(r+l,l;r2 f r+l) | / . |X(2,2;r2,2) H~ X(4,2;r2,4) ~t~ " " * ~l~ X(r,2;r2,r) • (A2) 

r 2 — 1 r 2 — 1 



We will now consider the case where the particles have 
n internal states rather than 2 internal states as they do 
for spin 1/2. Let us focus on the case n = 3, which would 
correspond to three-component states, such as the spin-1 
states considered in Refs. [H and[56| . The generalization 
to higher n will be clear after that. For n — 3, the 
starting product theory is Ad (3, r + 2) x M (r + 2, 2r + 
1) x A4(2r + 1, 3r). The bosonic field which one can take 
as the field which condenses is6=(3,l;3,3;l,3), which 
has scaling dimension hb = 2. As was the case for n = 2, 
we need that r mod 3^1, such that the factors in the 
product theory are well defined. 

Starting with r odd again, we can label the fields of the 
new (block-diagonal) modular invariant as follows. There 
arc a few different groups, namely (1, i\ 1, 1; j, 1), with i = 
1, 2 and j = 1, 2, . . . , 3r — 1, in addition to (1, i; 2, 1; j, 2), 
also with i = 1, 2 and j = 1, 2, . . . , 3r — 1. The fields 
appearing in the blocks labeled by (1, i; 1, 1) are 
(</i, i; <?2j 9i;i> 92)) where (71,92 mod 2 = 1. In the case 
of the blocks labeled by (1, i; 2, 1; j, 2), the fields appear- 
ing in these blocks are (q±, i; 92, qi; j, cfe); but now with 
<7i mod 2 = 1 and 92 mod 2 = 0. 

In the case that r is even, the labels of the blocks 
are slightly different, in the same way as was the case 
for r = 2. There are for different types of labels, 
namely (1, 1; 1, 1; j, 1), (1, 1; 2, 2), (2,2; 1,2; j, 1) and 
(2, 2; 2, 2; j, 2), all with j = 1, 2, . . . , 3r- 1. The fields ap- 
pearing in these blocks are of the form {qi, i\q2, qi', j, 92), 
where qi and (72 are either even or odd, depending on the 



type of block they belong to. In either case r even or 
odd, the number of fields in the block-diagonal modular 
invariant is given by 4(3r — 1). 

A few remarks about the case of general n. The prod- 
uct theory one starts with is M. (3, r + 2) x M. (r + 2, 2r + 
1) x ■ • • x M{3 + (n - l)(r - 1), 3 + n(r - 1)), where 
r mod 3/1. The boson with scaling dimension hb = 2 
has the labels (3, 1; 3, 3; . . . ; 3, 3; 1, 3). The number of 
different blocks, i.e. the number of fields in the new 
modular invariant, is 2™ _1 (2 + n(r — 1)). The precise 
form of the labels of these blocks, and the fields ap- 
pearing in them, depends on r being even or odd, in 
a way which should be clear from the case n = 3 above. 
The general form of the fields in the blocks is given by 
(Qx,i':Q2,qx',q3,Q2; ■ ■ ■ ; Qn-x, Qn-2',j, Qn-x), where i = 1,2, 
j = 1,2, . . . , (2 + n(r - 1)), and q 2 ,<73, ■ • - ,9n-i can ei- 
ther be even or odd. For r odd, q\ is always odd, while 
for r even we find that i and q\ have the same parity, 
(i + qi) mod 2 = 0. 

The invariants we discussed above are relevant for 
paired states with k = 2, based on particles with spin 
(n— l)/2, with the property that the wave function van- 
ishes as r powers when three particles with the same s z 
come together. It is possible to generalize this construc- 
tion to arbitrary clustered states with k > 2, but these 
will be based on cosets of the type su(k)\ x su{k)\/ su(k)2 
etc. making the structure of the modular invariants 
somewhat more involved. We will not deal with the case 
k > 2 here. 
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